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NOMENCLATURE 


V D = velocity of diffusion 
C = speed of light 
kT — temperature in eV 
p = kinetic pressure 

= electron temperatures (eV) 

= ion- temperature (eV) 
n = number of electron or ion-densities 
B = magnetic field in gauss 
n Q = neutral density 
a = plasma radius (cm) 
a> c = plasma frequency 
= resistivity 

L = length of plasma column 
= ion-ion mean free-path 

, g 3 

U g = average electron energy density <-10 eV/an 

V 3 

U ± = average ion energy density (10 eV/cm J ) 

9 

A q = obstacle area (cm ) 

P = microwave power (watts) 
f as microwave frequency (GHZ) 

Q q = unperturbed cavity Q q of plasma chamber 
t = time (sec) 

g 

thermal conduction time of species 
Tp = plasma column particle confinement time 
e s = average energy loss per species 
t__ = electron energy equilibration term 
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CHAPTER- 1 

The upsurge of our indutrial civilization in the 

first half of the twentieth century. was founded upon 

fossilfuels- coal and oil. But there are many "have-not" 

nation^ and even contries rich in this fuels are now seeing 

« 

rapid depletion of their reserves At this juncture uranium 
has come to the rescue as a hope for the future. The worlds 
uranium and thorium, it is estimated, represent an energy 
reservs,somewhere between ten and hundred times larger than 
coal. Even so fissinable fuels t<9©- are an exhaustible 
supply. At the rate at which the world's energy needs are 
expanding practically all the economically recoverable 
uranium could as coal, might be exhausted within another 
century 'or so. 

Besides limited uranium reserves, fission power 
also presents a more immediate problem; namely, disposal 
of its rad^active wastes. 

All this helps to explain the drive, if not race, 
to fiwid out whether thermonuclear power can be tapped and 
put to work. If the fusion feaction can be made to yield 
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power , it will solve forever both the fuel supply 
problem and the problem of radioactivewastes . The basic 
fusion element deuterium, is as inexhaustible as the 
oceans, and fusion produces no appreciable amount of 
radioactive by products. 

Nuclesr fusion is not exactly a new phenomenon . 

It has been generating the power of the Sun and other stars 
for billions of years. But to create and control fusion power 
on earth is a problem of a totally different order from 
harnessing fission. It is undoubtedly the most difficult 
project ever presented to scientists and engineers. 

The fission releases energy, because part of the 

mass of the fusing nuclei,.’ s transformed into energy 

2 

according to Einstein's famous equation E= me . A 
continuous fusion reaction, is analogous to the familiar 
process of combustion. To stick together or fuse, the molecules 
must collide violently, which means the material must be 
heated. Thrree condition are needed to burn a chemical fuel 
and harness its heat to do work; (1) The fuel must be 
raised to its ignition point. (2) There must be enough of 
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it to sustain a continous reaction. (3) The energy 
released must be tapped in a controled manner. Now 
precisely the same cord ition are needed to make a 
nuclear fusion reaction go and do useful work. The great 
difference is that for a fusion reaction the ignition 
point is rather high- hundreds of millions of degrees 
centigrade . 

This one condition - the attainment of which 
was quite unthinkable on the earth until recently - 
underlies all our problems. The. answer to the question 
"what is a controlled fusion reactor" is best given today by 
merely saying that it is a device within which appropriate 
isotopes of light elements could be caused to undergo 
nuclear fus^ion, the end result being the controlled 
production and extraction of useful quanitities of 
energy, in excess of that required to operate the device. 

Among the nuclear reactions which appear promising 
for a controlled fusion reaction are those involving the 
various isotopes of hydrogen, helium and lithium. Typical 
reaction are given Table (1). But unfortunately nuclear 



TABLE -l-j 


TYPICAL FUSION REACTIONS 


1. D+D 

wm 

3 He 

+n 

+ 

3.25 

MeV 

2. D+D 

= 

T 

+p 

4 

4 

MeV 

3. D+T 

as 

4 He 

+n 

4 

17.6 

Mev 

4 . 3 He+D 

= 

4 He 

+P 

4 

18.3 

Mev 

5 . 6 Li+D 

- 

a^He 


4 

24.4 

Mev 

6. ^Li+p 

= 

4 

a .He 


4- 

17.3 

MeV 
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fusion reactions have significant reaction cross-section 
at relative energies ranging between ten kev and hundred 
keV. At these temperature the matter is completely 
ionized and exists in the form of plasma in which the 
ion and electron densities are equal to maintain charge 
neutrality. 

There are basically two approaches to fusion. The 
first of these is magnetic confinement of plasma. The reacting 
plasma is too hot to be contained by mater ialw alls , and it is 
therefore proposed to use a "magnetic bottle" i.e a magnetic 
field configuration which curls up the particle orbits in 
such way as to greatly reduce their escape rate from the 
reaction volume. Speaking in a very general way, the confined 
state of plasma is far from equilibrium, and there are * 
many mechanismsleading to the break up, or dispersal, of the 
confined plasma. Research over many years has lead to an 
understanding of these escape reutes and many have been 
eliminated. There is , however one phenomenon which is 
not yet understood; it was discovered by Bohm in 1949, 
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According to him, there exists an amomalous diffusion 
process in a magnetic field B. The particles diffuse 
with a velocity given by 


V = - ckT ,—t 

D l6^n B V n e 

c 


8.64xlO^T < 7 - 

^B n e B 


where is a factor which theory cannot account for in a 
straightforward way and which according to Bohm is 16. 

Much research in recent years has consisted of attempts 
fco discover which circumstances are reponsible for this 
anomalous diffusions, presumably caused by turbulent 
fluctuation process. 

Putting, <Tn e = n^/r 

shows that Bohm diffusion limits the life time, T = r/v D , 
of a plasma confined in cylindrical geometry with a radius 


r , to a value approximately given be 

T = 10 4 (Y n V. r 2 B)/T 


( 1 ) 


where ~/j is an improvement factor which may be raised by 
1 a 

controlling the fluctuation level. There B is measured in 


.8 


Gauss and T in Kelvin so that for T = 10 and a field of 
hundred kG and r = 100 cm, T= 16ms , if V. = 1 . 
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Magnetic confinement systems are characterized 
by the ratio of thermal energy 2nkT and magnetic energy 


B 2 /8 , 


= 16 K n kT/B' 


The parameters of a possible reactor can be 

determined if we also consider the power loading of the 

2 

reactor walls. A figure of 50CW/cm is commonly taken as a 
reasonable value. For plasma contained within a radious r 


the number - N of particles per unit length is given using (1) 

and the power W reaching unit 


by N= r 2 n = 10 4 n / B 


area of a wall of radius R is ( using (2)) 

3NkT 3X10 14 B T 4. 

w = — Am = — ^ — 7m — / 0 = 


2 RT ~32 R~^ 


8 . 


Again for a temperature of T= 10 &,we find that the wall 


radius must be at least 


R= 5.9 fV B 


(i/-i 


(3) 


Now R must be larger than r and this condition 


yields/ 


B ^ 3 (“6 


1/5 


-£ /F Zi. kilo gauss 


(4) 
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For the minimum value of B when r is determined 
from above. 

A reactor with R=100cm is already a very large 

2 

device. Adopting this value and B = 2x 10 k G which is 

practicable, we see from (3) that must be at least 1300 to 

satisfy criteria laid in (3) . With this value and E =0.01 

P 

E must be larger than 180 kG which is consistent with our 
assumption of 200kG. Such values of have bot been achieved 
so far, even at temperatures smaller than thermonuclear 
onees but it may well be possible to reach this goal. 'We 
have presented these estimates to show that the search 
for alternative methods for obtaining fusion power is still 
necessary. 

The foregoing data imply confinement times of 
lms and densities of 10 /cm not easy to maintain for 
sufficiently long times for extracting power. 

This has led to an alternative scheme of confinement 
namely inertial confinement by lasers. It is for this reason, 
that microexplosions involving much higher densities during 
much smaller times are tried using lasers. We now tern to the 



8 


problem of heating and compressing a pellet containing 
a deuterium-tritum misture. We see it now as a spherical 
object with a radios of the order of 100 micrometers . The 
scheme is to illimunate this target as uniformcly as possible 
with intense laser light whereupon some light will be 
absorlced and materials will evaporate from the pellet 
surface. At high temperatures it will ionize and form a 
rapidly expanding plasma mantle. Under intense illumination, 
the plasma density may rise to such an extent as t o make it 
opaque to laser light. The heat from the absorption region 
will be conducted to the pellet where it will cause 
further evaporation. The pressure of the plasma so ' -gtflPmecl 
may become very high that a shock wave is launched into the 
target. The mass loss of the pellet by fast evaporation is 
called ablation. 

Thus accofding to our criteria, with E = 0.01, 

ir 

nT= 10^ , an explosion lasting lOOps would correspond 

O O 

to a density of 10^ /crrt which would be just sufficient 
provided that the fuel could be kept at this density 
during this time. 
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Following is a brief review of the various 
fusion devices being experimented on. 

THE PINCH Device 

The theory of the constructed gas current was 
first developed by W.H.Bennet in 1934. A few years later 
a different approach was presented independently by L.Tonks 
who used the term "pinch effect" in the sense in which 
it is employed at the present time. Subsequently , a number of 
other theoretical analyses have been published of the 
construction of a current filament due to the action of the 
azimuthal magnetic field generated by the current itself. 

■tv 2 2 

Using the Benet relation I = 2N g k (T e +T^)c it appears that 
the attainment of a sufficiently high temperature in a pinch 
discharge depends on the ablity to pass current of that 
magnitude. 

Since the temperature in the quasi-equalibr ium pinch 
varies as the square of the current* Hat equilibrium 
temperature is limited as the rate of energy loss will increase 
rapidly with the pinched discharge current. 
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A considerable amount of work has been done with 
"fast discharges" in which the constricted plasma is 
established in a time that is short in comparasion with 
the sound transit and the magnetic field diffusion times* 

The dynamic pinch ( or shock pinch) as described above 
has received a great deal of attention experimentally since 
it offers, in principle, a way to obtain high average particle 
energy in a plasma* When the pinch reaches its smallest radius 
the kinetic energy of the particles in the plasma reach their 
maximum inward velocity. If the pinched discharge column 
can be kept contracted long enough to permit randomisation, 
i.e., thermal ization of the ion velocity to take place ,a 
highc. temperature will result. 

One of the major probems associated with the 

constricted discharge and infact, with any system of 

charged particles confined by a magnetic field- is that of 

instability. There are two main types of instablity which 

leads to the destruction of a pinched plasma. The first 

often called the sausage type of instability, results from an 

azimuthally symmetrical localized constrAit-ion (and expansion) 
of the plasma column. 
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This is diagramatically represented in figurel-l 
Consequently, the constrictions tend to become more and more 
constricted , whereas the bulges will grow larger and larger 
until the pinched discharge is completely disrupted. 

The second type called jaink instability originates 
with the formation of a band or a kink in the plasma column 
while the latter retains its uniform , circular cross-section. 

This is d iagramatical-ty represented in figure 1 - 2 - 
It is seen that the lines of force of the azimuthal self- 
magnet icfield , due to the current in the plasma, are brought 
closer together on the inside, but they are further aP a rt on 
the out side of a bend. Thus once a slight kink developes ,it 
will grow in size until the pinched discharge strikes the walls 
and is cooled. The plasma will then become diffuse, and fills 
the containing tube. 

A part from the problem of instabilities, the theta 
pinch devices are also bogged by problems of end losses. This 
has le-d to the the development of toroidal devices. 
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The adiabatic in variance of the magnetic moment in a , 
space dependent mganetic field leads to a magnetic field * 

configuration known as a magnetic mirrors. 

The other concept led, to development of a toroidal 
pinch device initially. The most promising of all thermonuclear 
reactors to date is of the Tokamak series developed in the 
Soviet Union. This is a toroidal pinch with a very strong 
toroidal magnetic field B- generated by the current 1^ A poloidal 
field results from a large current induced in the plasma by an 
iron core. The stellarator was initially conceived as a toroidal 
like steady state confinement. The simplest form is that of a 
system of coils which generates a magnetic field in a toroidal 
region of space. By contradistinction with the toroidal pinch 
the magnetic line of force here are "parallels" of the 
tordida .instead of medians. In addition there are a number of 
other devices like the spheromak and base-ball machine. 

Recently , it has been found that in order to achieve a 
practical fusion device one has to use a high beta fusion 
device. Given the low-beta of tokamaks, there is renewed 
interest in theta pinch devices and more recently the 
REVERED -FIELD PINCH (RFP) . 
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The main plasma parameter of a thermonuclear 

configuration are determined by physical and economical 

consideration. For the output in Thermonuclear energy from 

fusion reactor to exceed the energy needed to heat the plasma 

the product of the plasma density and the time T should 

14 3 

satisfy the Lawson ceriteria; nT=10 2 /cm where T is 

the mean time a charged particle can be held in the plasma 
volume. The mean power density desired from a fusion reacter 
is determined by the condition that the capital cost of the 
reactor should be comparable with that of the conventional 
power sources. The maximum power density possible is limited 
by the capability of the srounding material to withstand the 
energy released, without being destroyed. Under these two 
condition, the power density released by the fusion reacter 
should be btween 1 and lOOW/cm . Xfc. is clear that power 
release will be a function of the temperature and the density 
of the confined plasma. Under the circumstances, it is possible 
to consider three regions of confinement: these are: 

(10 high density low temperature, short life tine 
plasmas. 
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(2) low density, low temperature, long life time 
plasmas . 

(3) low density, high temperature, lone life time 
plasmas . 

Whatever is the form of the magnetic confinement 
the plasma can always escape from the configuration. Energy 
losses from a controlled -thermonuclear reaction systems may 
be divided into two broad categories- direct losses and 
indirect losses. Among the direct losses are those due to 
plasma instabilities, radiation, charge exchange, thermal 
conduction and diffusion. These are all energy losses that 
take place directly from the plasma itself. The indirect 
losses are those arising from external circumstances, such 
as the losses resulting from joule heating in the magnetic 
field coils and from the conversion of the thermonuclear 
reaction energy into useful power. 

Some of the parameters of a thermonuclear reacting 
system, notably the dimensions of the reacting chamber and the 
strength of the confining magnetic field, can be yaried 
in such a manner as to increase the power 4 output and 



15 


minimise losses. The broad rules which are ^pplic able are 
refered to as scaling laws. . 

One of the main uses of a zero dimensional 
simulation program such as this might be to predict scaling 
laws for magnetic confinement devices. 


1.2 SOME USEFUL DS FI NT IONS; 


In this section, we first give definitions of system • 
parameters like beta, neutral pressure , crowbarring, elastic 
and inelastic collusion, ion an mean free path, and ambipolar 
diffusion. 


BETA;- it is often convelt&ent to determine the 
ratio of the kinetic pressure of the particle to the 
external magnetic pressure. The dimensionless quantity 


is then defined by 




B^/8 PC 


The energy balance equation is given by 


p4 _B^ 

8K 



8 



This may be written as 
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Since the minimum value of B is zero, the ratio 
has a maximum value of unity; this i^ould represent the 
ideal case of a perfectly diamagnetic plasma from which 
the magnetic field was completely excluded. In this event 


we have 


= f B r 


max 


8tT 


where p Q is the maximum kinetic pressure of a plasma that 

can be confined, in a study state, by a external magnetic 
A. 

field of strength B . 

NEUTRAL 'PRESSURE This is the pressure exerted 
by the neutral gas* in the system. The magnetic region in 
which the injected ion are trapped contain a large number of 
neutral gas atoms, owing to the finite initial pressure 
the system must have. 

CROW BARR INS ;•= This is a technique# by which 

the external magnetic field applied to the plasma is 

forcefully brought to zero after a predetermined amount of time. 

ION-ION MEAN FREE PATH:- This is the mean distance 

travelled by the ion between col&sions with oter ion. 
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PLASMA -CONFINEMENT TIME : This is the average 
time the plasma stays confined before it escape from the 
magnetic field. This is also equal to 1/e times the plasma 
decay time, 

AMBIPOLAR DIFFUSION: This is a diffusion process 
which takes place due to the movement of electron and ions 
across the magnetic field with equal velocities, 

1.3 LITERATURE REVIEW : 

Computer simulation of various fusion devices 
using time dependent zero dimensional computer models 
have been carried out by several authors (S,V D ). 

In reference ( 4 -) the author has indicated that 
the mod&^that he has developed was a fast /Simple/ time 
dependent , aero dimensional computer model for magnetically 
confine# plasmas. According to^author, this method could used 
to model other devices or heating methods very easily. 
According to the author most of the terms do not depend 
on the particular geometry, containing at most the total 
plasma volume , and so could be used unchanged in modelling 
other devices. 
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Several theoretical studies which describe and 
losses from theta pinch plasma column have been reported 
in the zero dimensional code of Green et al ( ' 5> ) 
the electron and ion temperatures were assumed equal and 
particle end losses was included resulting in a solution 
for temperature which was separated in spatial&time varialbles. 
The zero dimensional mod<2^ developed by Kievans et al 
includes both uniform plasma profiles as well as the effects 
of ion thermal conduction and finite thermal conductivity. 

In the present thesis work, attempt has been made 
to include the geometrical th<£ta pinch device as indicated 
in reference (2- ) into the model developed in refereace( ± ) 
to produce a zero-dimensional model of a pinch device 
where the preionzation has been done using a plasma gun* 

1.4 OUT LINE OF THE PRESENT WORK:- 

In chapter 2 a zero dimensional model is developed by 
defining the various energy loss terms. These energy loss 
terms have been space averaged in order to make than 
space independent. From this a set of equations, one each 
for plasma density electron energy density and ion energy 
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density have been defined. In addition in chapter 2 the 
energy terms of a theta pinch device have been derived. 

In chapter 3 we established the method by which set of 
equations are solved. We also determine .in chapter 4 by 
parameter adjustment, the optimal neutral pressure for 
obtaining the appropriate variation of temperature. As will 
be seen latter , any one of the following three situations 
can occur” in a pinch device: 

(I) Progressively decreasing electron temperature 
variation whereas a ? progressively increasing densities. 

(II) Progressively incereasing electron temperature 
variation whereas a progressively decreasing densities. 

(III) Gradually approached peaks of electron ion 
temperaurs and densities. 

These cases are obtained by proper adjustment of 
neutral pressures. Next the electron and ion temperature 
variation in case of adiabatic compression studied. 

Finally, the summary and the main conclusion of this 
work are presented in chapter 5 . 
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CHAPTER 2 


The computer program (SIMULT) calculates the time dependent, 
spitially averaged density, electron temperature and ion 
temperature in a cylindrical multipole device [ \ ]. The 
program is zero- dimensional in the sense that it predicts only 
spatially averaged quantities. It has been used mostly to study 
continuous wave microwave heated plasmas, but pulsed microwave 
and gun-injected plasmas can also be studied. 

The numerical method consists of solving three coupled, time 
dependent, non-linear, first order, ordinary differential 
equations, one for the p article density, one for the electron 
energy density and one for the ion energy density. 


dn _ s >n 

dt ~ 3 t ionisation 

(A) 


3 n 

9t diffusion 
(B) 


3n 

3 1 field decay. 
(D) 


9 n . 

3 t obstacles 
(C) 


( 1 ) 


dU 

6 

dt 


3 U 


3U 


3 U 


3 1 microwaves 3t excitation 3t ion collision 

(e) (f) (g) 

3 U 3 U 

e _ e 

3 1 bremsstrahlung ” 3t synchro ton 

(H) (I) 


dt 


3U 

e 

3t thermal conduction 
(J) 

3 U, 


( 2 ) 


9U i 


3t electron collisions 3t charge exchange 
(K) . (L) 

... 3 _Zi 

3 t thermal conduction 

(m) 


( 3 ) 



The first equation (l) describes the electron and 
ion densities in terms of ionization# diffusion# obstacle 
losses and field decay. The second equation (2) describes 
the electron energy density in terms of microwave heating# 
excitation energy# electron-ion collision; br emss trahliitt^ losses; 
synchrotron losses and thermal conduction losses. The third 
equation (3 ) describes the ion energy density in terms 
of ion-electron collisions# charge exchange and thermal 
conduction losses. 
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The computer results could not be compared with experimental 
values because all the required parameters in the pinch device 
were not available* 

The three equations are solved by the simplest method of 
successive iteration using the current value of each parameter. 
The electrons and ion-temperatures are incremented according to 

AT = (£U - T An ) /n 

0 e e 

AT ± = (£U i -T i An)/n 
An is change in n in At. 

The iteration step At can be adjusted to insure any desired' 
degree of accuracy. The initial conditions (n,T ,T^) can be 
set arbitrarily. 

In this study a time dependent magnetic field as specif '- ed 
below. 

-t/T -2 

B = B Q e ' sinwt 0) = 200 /sec x = 10 sec. 

In this study microwave power pulse«~specified in [ Jl ] 

A. 

was us ed . 

The neutral hydrogen gas is assumed to consist of thermal 

(T = T ) molecules. The initial value of the neutral density 
w 

is an experimental parameter which must be specified. 
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Reference [j,] indicates that the neutral (hydrogen) density 
is calculated at each time step from the pressure (p) as follows! 

n D = 3 22p exp[-3 .07xl0" 4 ■ -fy Q e " 17 ' /T e] 

T e 

+ 12.88p exp[— 1.7x10 ■■ "T/o e ^ ei] . 

r jy m f a 


The first term represents the thermal (0.025eV) neutrals and 
n FC re P r esents Franck Condon (7eV) neutrals. The exponential 
factors account for the finite penetration depth of neutrals 
into the plasma. The thermal neutrals are assumed to be at 
wall-tonperature (T t ) . The experiment is assumed to take place 
on a time scale short compared with the pumping time so that 
no neutrals are lost to tjSe vacuum pumps. When the neutral density 
is increasing in time, the rate of increase of each component 
is limited to 


dn TH ^ 1.9xlO ~ 5 


dt 


*n 


TH 


^FC > 
dt 


3.4xlO~ 6 n 

a 


EC 


representing the transit time of a neutral from the edge to the 
center of the plasma. 

(a) IONIZATION -The ionization rate for a Maxwellian electron 
distribution in a cold neutral gas, can be written in a form 
given by Drawn [ref(£ )]• For hydrogen the electron neutral 
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particles migrate across B by a random walk process. When an 
electron collides with a neutral atom, the electron leaves the 
collision in a different direction but its phase is changed 
discontinously. 


In such a case, D . = =— 

1+6) T 


When u 2 T 2 >> 1, D ± = 


KT 

inv 


KTy 


2 2 
03 T 
C 1 


mo3 


63 


C 


=> D 


eB_ 

m 

2 

KTxn xcrv x m 
o 

2 2 
mxe xB 


KTxn o 

B 2 xe 


2 


<5v ysxi 


T 


j 


- lit x n 

-KTxn 0 x<7vxm 
= B 2 xe 2 


Kxnn T 
o e 

n 2 2 

B xa 


x 


8 n 
8 x 


Consider the I term • 

V xnxSkT T/nT 

In a fully ionized plasma, D 2 = “T“[ T eV ] 


B 


B 


T) xn n 

Electron current t = 5 x T 


KxT) , mx n 
n i 

B 2 

Since = constant times (kT) 


- 3/2 V&Lv&al 1 


" 0 



24 


collisional ionization rate is 

+ In (1.25 (1+1/S) ) ) . 
where S = 15.6/T . 


dn 

dt ionization 


- 371 e ~ S ( 
1+S T 

e 


1_ 

20+S 


From another approach the ionization rate for a Maxwellian 


electron distribution is approximated by an anlytical fit 
for the data given in [n]. In view of the inavailibility of 
ref[ n ] a analytic fit of the cross-section for ref [ l2 -] was 
performed and a similar variation was found. In general the 
ionization is largely proportional to the electron temperatures. 


nn 


dn _ cn . -“o -17T 

3% ionization * ^,1/8 6 e 


(B) DIFFUSION- Classical diffusion due to electron-ion and 
electron neutral collisions is given by 

2 0 nn T 

dn _ n n ir -3 o e 

dt diffusion ~ 2 2 m l/2 ~ xu n 2 2 

Dai' d a 

e 

(i) cn) 

The above term has been satisfactorily derived as follows’. 

Consider the II term. In a magnetic field charged particles 
will move along B as 

tj = nVj = ■+ M-jnE :-D7f n . 

If there were no collisions, particles would not diffuse 
at all in the perpendicular direction. When there are collisions 
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-3/2 

- „ constant x nxn_xT - xT 

3 n _ o 

at diff " «2 2 

3 xa 

constant xnxn 
__ . o 

B 2 xa 2 xT 1 ^ 

Here ambipolar electric field has been considered* 

(c) OBSTACLES: Hoop supports and probes which intercept 
the plasma are assumed to collect particles at a rate given 
by the ion saturation current for a isotropic plasma. 

Current density J = nx velocity 

= constantx ^(T^+T^) 

2 

A /a is a factor which determines the fraction of the 
area occupied by the obstacles. 

3 n 2.0xl0 5 nA c V CV^O 

Hence we have -g-g obstacles - T : 

a b 

(D) FIELD DECAY : When ^ is negative field lines are 
leaving the machine, and plasma is assumed to leave also at a 
rate given [ | ] 


3 n | 

3 t field decay 

(e) microwave heating : 


-min[0, |.B ] 

Ref[ \ ] calculate the microwave 


heating rate as 


31 " 

— §| 

3 t 1 microwaves 


2xl0 9 R 
(1+nyni a 2 L 
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where n c is the density above which the microwaves are 

totally absorbed. In general/ according to ref [ \ ] n is 

2 

given very roughly by n c = 487f /Q q . 

Here is a brief discussion about the same. The incidence 
of a microwave on a plasma gives rise to the following phenomena 
(l) Absorption (2) Diffusion/ reflection and secondary emission 

(3) Phase lag, caused by differences in refractive indices. 

( 4 ) Modification of the polarization and birefringence among 
other phenomena. 


Here Q- factor is the measure of the sharpness of resonance. 


Instead of a quantitative derivation of the expression on 
the R.H.S., a dimensional analysis in order to understand the 
derivation of the various terms is given below. 

Consider R.H.S. W e have - ene £2I - L.H.S. 

cmx( factor) tim exarea 


A little more meaning could be given if one considers the 

process to be due to transit time magnetic pumping. If it is 

assumed that there is equipartition of energy between ions and 

d.T 

electrons/ P = 3nk assumes that the phase velocity 

G)/k for optimal heating is of the order of the mean ionthermal 
velocity. 


“V 


Here K = 


2kT. 

<ir> 


P = 3nk 


dT 

dt 


P = K f 

dt 


( 1+n^/n ) a L 
2x10 
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(F ) EXCITATION An electron can also experience a large 
number of different types of inelastic collinsions with atom, 
ions and molecules. In a collision with one of these heavy 
particles, an electron may merely excite the heavy particle. 

The cross-section will be dependent on the energy of relative 
motion and the particular quantum level excited. Reference[ , X ] 
has assumed that there is a 30eV/collision. Hence we have 


9U 

— _ _ oQ S2. t 

3 t excitation dt ionization 


(G) ION COLLISIONS . Consider ^encounter between two 


charged particles. ■•ere each particle moves in a hyperbola 
relative to the m of the two particles. The distance of 
closest approach is p ''the impact Parameter''. Consider a 
test particle with mass m, a charge Ze/c and a velocity W . 

Then AW represents the charge of velocity of the test particle. 


Consider the Maxwell -Boltzmann distribution 
f (o) (w) = _nLi_ e -L 2 M 2 

where n and m are the particle density and the mass of the 
particle in question, and L is defined as L 2 = ~~ - 

To measure the rate of diffusion in the direction, 
we set N equal to the average number of encounters in one 
second. The resultant value of (AW^.) , measuring the 

increase of velocity dispersion of a group of particles/second 
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may then be denoted by < (Aw^) ^ . This quantity is known as 
7 ' diffusion co efficient' * . 


When the distribution function of the field particles 
is assumed to obey a Maxwell -Boltzmann distribution only 
three independent diffusions need to be considered, < AW^ > . 

^ ^ W ll ^ • > / an< 3- < (^W )j__ > . The first quantity represents 
"the rate at which moving test particles are slowed down by 
interaction with the field particles. The quantity < - (AW^)^ 
represent the rate of increase of (AW)^ in a direction parallel 
to the origional motion of the particles. The corresponding 
quantity < (AW ) l _ > represents the rate of increase in 
the perpendicular direction 


h/p o where h is the Debye Shielding distance and 
p Q the impact parameter 

3 


2ZZ^e 


.(££_> v 2 

iin ' 


As In increases T increases,* InAiecr eases as n increases. 


Detailed computation of the diffusion coefficients have been 
carried out by Chandrasekhar [ 3 ]* The resultant formulas 
for the three diffusion coefficients are 

<AW lX > = -A d l J( 1+ j~)G(L,W) 

9 A-. 

<(AW 11 )^>*s jp6(L ( W) 

A 

< (aw ) 2 > = ~ (l,w)-g(l,w) 
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where the "diffusion constant" ^ is defined by 




4 0 0 

8 n e njZ^Z^ln 


m 


^(x) is the usual error function, 

* ( * } = Oi7T~ 

A ' o 

and the function G(x) is defined in terms of (x) by the 
relationship 


G(x) = - £l (x) 

2x 2 

The "time of relaxation" denotes the time in which 
collisions produce a large alteration in some original velocity 
distribution. If one considers certain times which are defined 
in terms of the diffusion coefficients, we have 7 t^ 7 , 7 t^ 7 
is the time between collision in which reflection gradually 
deflect the test particles by 90°. i n other words 

< (AW ) 2 > tjj = W 2 . 

An energy exchange term t £ may also be defined by the relation 

< (AE) 2 > = E 2 

E 

The change of energy, AE in a single encounter is given by 
AE = | { 2WAW 1 +(AW 11 ) 2 +(AW ) 2 } 

Spitzer has determined the rate at which equipartition of 
energy is established between two groups of particles having 
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The actual calculation is difficult because V and V are not 
well known. Hence we as same that on electron of speed v has 
an impact parameter b with respect to an ion of charge 
which is assumed to be infinitely massive. The electron 
acceleration is 


2 q e q i 


4 TC £ mb 
o 


2 2 

with v << C and neglecting the cross-product/ we have 


dU 

dt 


2 * 2 
q v 

6 n £ CT 
O 


dU 

dt 


2 4 
q. a 
~e 


96TT 3 £ 3 c 3 m 2 b 4 
o e 


The total energy loss rate by the electron colliding with ions 
•at all impact parameters is 


, T q?q 4 (n . v ) 

dU tL^e i e 


dt 


max 

/ 


96 71 3 e 3 C 3 m 2 b . b' 
o e mm 


%2?fb) C 


Taking b max = » and b mln = 


2 n mv 
e e 


2 2 

q . q n. n 
M r^e r e 

x 24 JT e 3 C 3 m h 


8kT 


w = 


m 


= 48xlO“ 37 Z 2 n i n e T-V2 


Watt/cm' 
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Maxwellian velocity distributions, but different kinetic 
temperature. T and . Next he used (a) and 'averages this 
over a Maxwellian velocity distribution. We have 


dT 

dt 


T x -T 

t eq 


where t is the time of equipartition given by 

“Si 


3mm-^k 


3/2 


eq 8(2H) 1 ' /2 n 1 Z 2 Z, e 4 ln m m 


( ?+ .^) 3 / 2 
1 


This implies 


dT (T e - Ti )x 8(2n) 1 /2 niZ 2 z 2 e 4 ln 


Now 


dt 


3/2 , T ^ gn3/2 


3mm, k' J/ ^ ± 


m m n 


fa , n dr 

dt “ 3 dt 


dU 

6 

dt 


Kxn 2 x(T -T. ) 


e l 


T . -e-,3/2 


( ~ + ~) 
m m 


dU 


which compares well with = 


61n (T -T. ) 
e l 

*2/3 


(h) BR EM S S TR AHLUNG ^ According to classical electromagnetic 
theory, the total power U radiated from a charge q moving - 
with a velocity V and acceleration v is 


dU 

dt 




67T E C 3 (l-V 2 /e 2 )' 







(l) SYNCHRO IRON RADIATION t Again according to basic 


electromagnetic theory. 


dU _g; 
dt 


6 71 £ C' 
o 


_ 2 
In this case (v/C) << 1 


v 2 -(v x V) 2 /C 2 
(l-v 2 /C 2 ) 

and v = . 


The electrons are slightly relativistic. Their angular 
frequency of rotation w is determined by their relativistic 
mass, so that 


u o = 


q B[l- (v 2 /C 2 ) j 1 / 2 


m 


■ i / -| 2 , _,2 \ 1 /2 

=1 ll-v /C ) ' 


Subs tiu ting v^ = co o v we get 


dU 

dt 


q 4 v 2 B 2 
e 


6 n e m 2 C 3 [ 1- (v 2 / C 2 ) ] 


Expressed in terms of momentum. 


dU _ 
dt 


4 32 


6 it e m 4 C 3 
o 


/ 


which is to be integrated over the relativistic distribution 
f(p). Since f(p) is isotropic and since these are two direction 
associated with p but only one associated with p 1 ^ , we 
have 


2 


-2 
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The exact answer appears in terms of the modified Bessel's 
function (z) of the third order* The integral definition 
of the Hansel functions (z) of first kind and second 

order and the derviative relations for these functions 
are used to evaluate these integrals 

,4tj 2 n kT K (mc 2 /kT ) 

w = Cfi-J -3—- §- 

3 71 e m c me K 0 (mc /kT ) 

O 2 ' 


From the asymptotic expansion of the modified Bessel's 
function. 


2 . 2 
e 0)0 

W = » 

c 3 71 £ c 


n kT 
/ c e 

. v, _ 


_ kT 

■) (1+ | s* 


= 6..2xlO“ 17 B 2 n T [l+(T /204 ) ] 

GO 0 


Watts / cm 


Here 0)^ is expressed in terms of B, 

(j) THERMAL CONDUCTION 1 Heat conduction due to classical 
collisions and to obstacle losses [ y] is given by 


dU a „ 

§8 _ 2.5T 

3 t thermal conduction e 3 t diffusion 


+ 5T-54T 


3 n 

e 9 t obstacle 


(K) et.eCTRON COLLISIONS I Ions are heated by classical 
collisions with electrons. 

a"! au e 

dt electron collisions ” dt ion collisions 
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CHARG E EXCHANG E ! Reference [ i. ] determined an 
approximate analytic form for change exchange losses , obtained 
by averaging cross-section found in ref. [ i2-] 

3 U. 0.0136nn T 2 

i. _ o 1 

at ~ t ± +ioo 

Our first study had been to determine the exact expression 
itself. Later a sixth order polynomial was fitted for the 
data in ref [ o] and found similar variation in the cross 
section . 


(m) ION THERMAL CONDUCTION l Heat conduction due to 
collisions and to obstacle losses [ 7 ] is given by 


3U. 

i I 

3 1 thermal conduction 


2 *^ T i S' t 'diffusion" 1 " 2 T i ^ t 'obstacles 


Che first equation (l) describes the ^l^ctron and ion 


is of ionization, diffusion, obstac 


is oj 


densities in 
and field decay 
electron energy density in 
excitation energy, electron^ibn co 
losses, synchrotron l&^ses and thermal cb 
The third ' aqua 
in terms 

cond«ction losses. 


5n 


losses 

p_econ d equ a t i on^(2^^esdr ibes the 

xcrowave heating, 
j^sion, bremsstrahlung 
ction losses, 

(3) describes the ion en ergy^dehsity 


ion- electron collisions, charge exchange and>thermal 


The above formulation of ref. [4_] was modified in 2.2 
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to include plasma contraction terms as proposed by ref.[Z--]. 
The following terms (N) and ( P) contribute to the electron 
energy density term, whereas the term (0) contributes to 
the ion energy density term. 

2.2 ADDITIONAL TERMS FOR A TH5TA - PINCH 1 

(N) A certain amount of work is done on the plasma by 

the system when the plasma contracts contributing to the 

increase in electron and ion temperatures. The following 

dV dA 

derivation is done for the term P mr = nT L(-r— 2) . 

s dt s dt 

For the theta-pinch device described in Chapter 1 we 
have the radial pressure balance as 


B 2 (r ) = B?(r)+ 4 .03 xlO -2 °n (r ) (T +T, ) 
ex ex 

Here we have assumed that the radial inertial effects 

can be neglected and that the magnetic field has a small 

B 

curvature. Let £ a • 

a ~ 

e 

Assuming a uniform magnetic field profile with 



n(r ) = 


( 2.48xlO J -^[BV(T +T. )](l+£j) r< a 

1 e e i u _Q 

l 0 r> a 


According to ref. [ tA n = nCr) for a uniform profile. 
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3 2 

n = 2.48 x10 19 ( ^S_™ t ( 1 ^2 ) 

3<P 71 3B 2 I 

Using the expression = 2na: _ 

a m- _ o r r =a 

o 


where 


<P = 1 B . (r ) dr). 

A 

P 


Combining equation [ ] and [ ( 3 )]/ the desired expression 

as 


for 


dt 


is obtained 


as. 


dt 


i aA , dB- 

_ Y (2 — — + 2 — §) 

b a^A_ at B dt 

P e 


From another approach it can be established that 

< 2 a 1 -S a „ dB dA T 

a. 6^/ **, 2 0.51 P. 1 r rr\ / 1 


at ~ 2$. 


2a where ct— — — S 4. §. 2 — . — _R ^ 2 — r (_± ) 

,mjre B dt H 3 A dt T +T. L l e K 2 ' 

e P e 1 x th 

+ T i (b-)h(0.25L-A 11 ) + §■ 3- tE e +e id 


•th 


'P 


From the above two equations, for a uniform profile, we 
have 


QA I 

2 _ __ : 

A p dt 




dB (1-S B ) 


4* 


, CA y 2\ V B dt T +T. 
( 5 +S ) e e 1 

3. 


-A) 


where A = 


% — (e +e.)+T Q (~ — )+T. (-5 — )h(0.25L-\. .) 
3 t o 1 e a 11 11 


' th 


'th 


(0) A tefms similar to (N) above is used in the ion energy 
density expression* 
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(P) OHMIC HEATING : 


The ohmic heating term ft = L / rj j 2 (r)dA. 

y A 

A P 


With the assumption that the plasma is fully ionized the 
plasma resistivity 

T? = 3.27xlO" 9 ln /T 3 ' /2 [Ref (^) ] 

J- ^ 

9 B(r) 


j (r) from Ampere's law is I .7 

*0 ar 

For a uniform distribution = J — [B (0.5+0.5 (-) 2 ) 1 

O 2T d 37 ^ e 3. J 

0.5B x2r B xr 
_ • _ e _ e 

~ ™ ' 2 “2 
a a 

r _2 2 

p> B x~r 

fi e = L / 3 .27xlO" 9 x2.2/T 3 ' /2 x x 2 71 x r x dr 

O 3. 


Here 


r = ' r hZ/2ft . 

P P 
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CHAPTER- 3 

PROGRAM DESCRIPTION 


3.1 The numerical method consists of solving three 
coupled, time dependent, non-linear, first order, ordinary 
differential equations, one for the particle density, one 
for the electron energy density and one for the ion energy 
density. 


The three equation are solved by the simplest method 
of succesive iteration using the current value of each parameter. 
The electron and ion temperatures are incremented ac«ording 
to; 


T e = (A. U e -T e A n )/n 


T ± = (AU^MiVn 


The iteration step t can be adjusted insure any 

desired accuracey. The inita.l condition (n,T .T.) 

0 x 

are given as input and can be set arbitrarily. 


3.2 The calculation is done in the following manner. 

Intially we specify the peak microwave PO in wats , UMAX 
the number of iterations, DEA the initial neutral density. 



TABLE 3- i 

LIST OF DEVICES SIMULATED 


BIGLEV BIGUP 


TMAX 

0.043 

0.043 

A (cm^ ) 

50.0 

50.0 

AL 

800.0 

800.0 

AO 

o 

♦ 

o 

700.0 

QO 

20000.0 

20000.0 


LITTLE 

QUADX 

PINCH 

0.005 

0.002 

0.005 

18.0 

6.0 

18.0 

270.0 

160.0 

800 

90.0 

3.0 

90.0 

2000 

500 

2000 
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TEA the initial electron temperature , TIA the initial 
ion temperature and TMAX the duration of the experiment in 
seconds. We also specify the plasma radi^us , length of the 
plasma, the obstacle area, magnetic field B6, microwave 
cavity Q and the microwave frequencey F in GHZ. We specify 
the time step DT, the maximum neutral density DEMAX, the 
maximum electron temperature TEMAX, the maximum ion temperature 
TIMAX, thei minimum value of neutral density DNM, the maximum 
value of the magnetic field the microwave power , the ion 
saturation current and the maximum value of energy loss. 

We calculate the magnetic field and its derivative 
using the appropriate formula. The energy loss , BETA , the 
change in the neutral density DNT, are also calculated. After 
this the above calculations the current values of time DEA, 

TIA, B, P, SJTI,ELOSS £ BETA, are stored. 

Next the densities of electrons, ions and neutral 
atoms are incremented using the first order method. Following 
this wq '.check for the time taken and the program is allowed 
to run or terminate appropriately. 



FLOW 


CHART 
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3.3 SYSTEM PARAMETER ADJUSTMENT 

MAGNETIC FIELD:- Any time dependent magnetic 
field can be specified. In this study we use 

B(t ) = sin oit 

MICROWAVE POWER:- Any time dependent microwave 
power pulse can be specified. Reference ( i, ) has used, 

P = P Q MAX (O, 2B/B -1) 

Though the initial conditions could be set quite 
arbitrarily/ . . in this study the initial conditions have 

been set according to references ( J, ) . A very interesting 
outcome of the simulation is regarding the independence of 
the plasma parameter variation with initial electron and 
ion temperatures on the one hand and their sensitivity the 
neutral pressures on the other. 

3.4 ADDITIONAL PROGRAM FEATURES:- An interactive graphics 

package has been incorporated in the program to obyiate the 

time consuming excercise of plotting the parameters on paper. 

In addition to the large, amounts of time saved, it also 

ensures greater accuracy in flotting. 

The program is adequately representated in the flow 
chart given in Figure A) • 
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CHAPTER-4 

4.1 This study simulates a number of multi pole devices and a 
theta-pinch device whose dimensions are given in table [3-1]. 

The first four devices correspond to multipole machines 
simulated by the Wisconsin group [ i j . The pinch devices 
is similar to SCYLAC [ 2-] • 

The intent of this chapter is to simulate the time 
variation of various plasma parameters like the temperature/ 
density and Beta. With the help of the zero-di/ ensional model 
of the magnetic confinement system obtained in Chapter 2, 
a detailed dynamic analysis of the system for various neutral 
pressures has been carried out. The main theme of this 
process is the prediction of the plasma scaling laws and 
determination of the parameter space in which the system 
may be made to operate. It is important here to recall the 
fundamental assumption associated with the zero-dimensional 
model, that all the variables in the computation can be space- 
averaged over the cross-section. The second assumption is that 
the resistivity of the plasma can be computed with the assumpt- 
ion that the plasma is fully ionized. This can be somewhat 
justified as one finds that the density of neutral atoms-- 
remains at a constant initial value though the ion and electron 
densities grow to 10-15 times the neutral density. The third 
assumption is that the plasma confinement time T p is a variable 
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dependent upon the interaction of neutral atoms with ions 
and electrons. The question of how correct are the above 
assumption is difficult to answer with any generality, and one 
is usually reduced to verifying the validity of the assumpt- . 
ions in various specific^- The sacrifiosin curred in the 
process thus limits its validity to a certain range of operating 
conditions. The most important aspect of these numerical 
studies is to determine whether all the experimental observations 
can be verified using the existing theory. For this we analyse 
the zero dimensional model by numercally solving the set of 
linear equations for predicting the transioit behaviour. 

This qualitative comparision helps us in confirming the various 
types o f processes which are assumed to take place. 

4.2 THETA- PIN CH l 

Fig .(4*2 ) is a plot of the time variation of the electron 
ion densities. The time variation of the electron or ion 
densities indicates that higher values of densities of electrons 
and ions are obtained at higher neutral pressures. This is due 
to large contributions of the ionization term. One also notices 
that the peak keeps on broadening, indicating a shift towards 
temperature equilibrium. There is also a noticeable shift 
in the peak or a delayed peaking with decreased neutral density 
indicating that the ionization term becomes effective only 
towards the end with the decrease of the loss terms. 
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Fig. {4.3) is the time variation of electron temperatures. 

The time variation of electron temperatures indiciates the 
strikingly large values of temperatures that can be obtained 
using the pinch effect. The temperature variation also 
indicates a peak broadening indicating an early temperature 
equalization. This can be attributed the smaller number of 
neutral atoms present initially. The shifting of the occurance 
of the peak towards the starting point indicates the predominance 
of electron energy loss terms with increasing neutral pressures. 

Fig* (4.4) Is the plot of time variation of ion temperatures. 
One observation which one can make on seeing the time variation 
of ion temperatures is that the was do not gain significant 
energy (_ 60e”) . This is primarily due to the fact that the 
electrons being more mobile, absorb all the in coming microwave 
power and there is little time for ions to gain energy via elect- 
ron ion collisions. One also notices that the maximum ion 
temperatures occur at neutral pressures which are neither too 
high nor too low. Rather they occur around temperatures 
7 which are 7 optimum 7 as will be seen later on. 

Fig. ( 4 . 5 ) is the plot of the largest peak of 'BETA' in 
each of the simulation runs. Contrary to usual expectations, the 
BETA value goes through a minimum somewhere around the 
'optimum 7 values of neutral pressures. The BETA values tend to 

be very large either for low neutral pressures or for high 
neutral pressures. 
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Fig. (4.6) is a plot the variation of the final plasma 
cross-section with various neutral pressure. The steep drop 
could be due to neumerical instability. 

Fig. (4.7) is a plot of the plasma cross-sectional area 
as a variation with the operating time. The gradual increase 
in the cross-section is due to increased electron and ion 
densities and temperature resulting in decreased mean free path. 
The decrease in the cross-sectional area is due to decreased 
electron energies towards the latter half of the experiment. 

The final cross-sectional area is very law in sharp contrast 
to other values in the experiment primarily due to very law 
electron densities and temperatures. 


Upon gradually varying the values of neutral pressures 

12 

such that the initial neutral densities vary from 1.6 x 10 
11 3 

to 6 x 10 particles/cm it was found that one could obtain 

reasonable variations of densities and temperatures only for 

neutral densities varying from 5.1x10'*'^ to 5~. 8x10^ particles/cm^. 

Upon assuming neutral densities either below 5.1x10"^ or above 

5.8x10 particles/cm the density variation or the electron 

temperature variation become insignif icant. In this study 

11 3 

neutral densities of 531x10 particles/ cm was used as shown 

in fig, (4.2)/ (4.3)/ (4.4). 
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Fig. ( 4 . 8 ) is the plot of normalized values of plasma 
parameters with time at optimum neutral densities (531x10"^). 

One notices an early and broad electron temperature peak. 

This is because the microwave power input is primarily absorbed 
by the electrons initially. As predictable the magnetic field 
variation is a sinusoid peaking somewhere midway. The electron 
or ion densities shows a peaking in the latter help of the 
simultation. The ion temperatures show a peaking midway 
between the electron and density peaks . There is only a small 
space of time in which the ion temperatures, electron temperat- 
ures and densities show a peak. Significantly though the elect- 
ron and density variations show a broad peak, the ion 
temperatures shows a sharp peak. The time lag of the B peak and 
density peak is noticable in the Wisconsin work ( d. ) . 

OCTOf’O L6. 

Fig, ( 4 . 9 ) is a plot of electron or ion density variation 
with time for various neutral densities. Unlike the theta- 
pinch , the density peaking occurs at the same location and 
there is apparently no shifting in the time of peaking. 

Secondly the peaking is also not as broad as that of the 
theta pinch. This is because the ionization tern becomes 
predominant only midway in the simulation. 

Fig. (4.1Q> is a plot of electron temperature variation 
with time for various neutral densities. One notices the 
increased maximum values of electron temperatures with 
decreased neutral pressures. There is also a predominant 
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tail coruing into the picture at low neutral densities. This 

giv:2S an indication that the losses which were predominantly 

* 

due to Brammstrahlung and synchrotron were significant only 
then. There is also a considerable amount of peak broadening 
indicating that ideally the electron temperatures attains a p$ak 
value and retains it till the end of the simulation. 

Fig. (4.11) is a plot of ion temperatures varying vgfflTfftog 
with time. The ion gain temperature only in the latter half 
of the simulation indicating better electron-ion coupling then. 
In fact in the last f ew time steps the ion temperatures show 
an almost asymptotic increase in values. Even then the ions 
attain values of around (80aV) at very high neutral densities. 

Fig. (4J2) is a plot of the variation of BETA with time 

at various initial neutral densities. It can be observed that 

maximum value of the BETA peak occurs at neutral densities 

9 3 

which are around 17000x10 particles/cm . This slight deviation 
is due to the fact that (nT) product maximum occurs at this 
neutral density. 

Similar to the theta-pinch, it was determined that neutral 

12 1 

densities around 1.7x10 particles/cm were most suitable. 

Fig. (4,14) is a plot of normalized plasma parameters with 
time. . Unlike the case of a theta pinch fig. ( 4 * 8 ) the electron 
temperatures and densities not have a broad peak. As one would 
have also noticed the maximum values of the various parameters 
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are several orders less than th@se corresponding^ the theta 
pinch. Also present in the ion and electron temperature 
variation is the tail which attains very large values towards 
the end of the experiment. This is primarily due to reduced 
losses towards the end of the experiment. 

Ref. ( ^ ) has indicated that the energy loss rate can be 
taken to be proportional to -dB/dt. Upon incorporating such 
term one can obtain the plot of energy loss rate varsus time 
as shown in fig. ( 4 J3) . Ref. ( i- ) has indiciated that this 
can be taken to be the rate at which the magnetic field lines 
leave the system. Fig. (413) indicates that the energy loss 
rate grows rapidly with time reaching very large values towards 
the end of the experiment. There is no energy loss in the first 
half of the experiment because here the magnetic field lines 
grow stronger with time* 

THEORETICAL BASIS FOR NEUTRAL PRESSURE VARIATION t 

In this study the pinch formed after the pre-ionization. 
Although heating is not a problem for .the high energy injection 
device, a new difficulty arises. The magnetic region in 
which the injected ions are trapped contains a large number 
of neutral gas atoms, owing to the finite initial pressure 
which any systan must have. The neutral gas produces a steady 
loss of the trapped ions by the process of charge exchange. 

In this process the hot ion captures an electron from the 
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gas atom and thus escapes from the system as an energetic 
neutral atom. Thus charge exchange competes quite seriously 
with the desired nuclear reactions in the gas. 

It would thus seem that the first major barrier which must 
be overcome by a high energy injection device is the neutral 
gas background. 

DIRECT (BRUTE FORCE) PUMPING ! 

The most obvious method, of course, is simply to reduce 
the initial pressure of the system until charge- exchange no 
longer competes with the usual plasma processes of Coulomb 
scattering and nuclear reactions. The loss rate of ions/unit 
volume can be shown to be limited by the charge exchange. 

NEUTRAL ' * BURNOUT 7 ' * 

It is possible to obtain a condition for creation of a 
thermonuclear plasma which is less restrictive than the one 
above. This takes advantage of the fact that the self 
shielding action of the plasma tends to lower the neutral 
atom density in the plasma interior. Neutral atoms are removed 
in two ways. One is by the charge- exchange process itself 
since the fast neutral which leaves the plasma region. The 
second and' more effective, process is simple ionization of 
the neutrals by the fast ions (and electrons) of the plasma* 

As long as the rate of destruction, of neutrals in the 
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plasma interior is small than the rate of influx of 
neutrals from external vacuum region, there will be no 
appreciable difference between the neutral density internal to 
the plasma and external to the same* There is a critical 
point, however for the injected current at which point the 
neutrals are being destroyed as soon as they enter the 
plasma region. Once this point is passed, the neutral 
density drops in the interior and the plasma density increases, 
hence burningout more * 'neutrals' ' . This critical point 
will be called ''burnout'' . It is expected that the neutral 
density will drop sharply from its vacuum value when the 
injected current is only slightly above the burnout point. 

A crude first estimate of the burnout current is 
obtainable directly from the considerations above. 

TRAPPED ION PUMPING (GUTTERING) : 

A certain fraction of the neutral atoms of the system 
become ionized upon passage through the trapped plasma region. 
These ions are then constrained to move along a field line 
and will drift out to the ends of the mirror container. If 
suitable pumping arrangements are then provided, this 
mechanism could represent an important addition to the total 
pumping speed. In some cases it may dominate and appreciably 
lower the entire pressure in the interior of the system. In 
this case the trapped ions, of course, are acting exactly 
as an ion-pump, or as a getter. 
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It has thus been demonstrated that several aspects of 
the dynamics of pinch formations can be brought out using 
even a zero dimensional model of the pinch device'. 
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CHAPTER-5 

5.1 SUMMARY AND CONCLUSIONS: 

In the present work we have developed a zero-dimensional 
model of a octopole and a theta-pinch device for simulating the 
dynamic behaviour and properties of various plasma 
parameters . 

Chapter 1 is devoted to a general description and 
some useful definitions related to the present work. It also 
contains literature review and outline of the present work. 

In Chapter 2 we obtained zero-dimensional models of a 
theta-pinch device and a octopole by deriving the various 
energy loss or gain terms. Our aim was to study the time 
behaviour of various plasma parmaters under various initial 
operating conditions, such that one could determine if 
possible parameter spaces of feasible and stable operation. 

A very acceptable variation plasma densities and 

temperature is is obtained at neutral densities between 

9 9 3 

580 x 10 and 510 x 10 particles / cm for a theta pinch. 

9 3 

If the neutral densities are higher than 580 x 10 particlfe/cm the 
variation of temperature s becomes very predominant with 



the electron and ion densities tending to be very low. 
The variation of the plasma cross-sectional area is also 
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very typical of a theta-pinch device. The variation in 
temperature and densities for the octopole is several 
orders less compared to a theta-pinch device for a given 
set of initial conditions. 

5.2 FURTHER SUGGESTED WORK 

The first improvement that can be done cancer ns 
the method adopted to solve the equations of the zeroT 
dimensional model. In this study a first order method has 
been used. More accurate results can be obtained if the 
equations were solved using a second order method (Appendix a). 
J.R.McCowan et al (I 5 ’) have indicated in a paper that 
Heindrich et al (^) found the rotational energy to 
be comparable to that of the thermal energy. Hence the 
code could be modified in order to be able to examine 
possible effects of rotation on transport phenomenanin 
the plasma. The zero D model includes five first order 
differential equations. Four of the equation are derived 
from expressions for the time rate of chance of electron and 
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ion temperatures, total particle in ventory and the 
magnetic flux within the*, r plasma. The fifth equation is 
obtained from radial pressure balance. 

Another extension could be in the direction of 
developing a zero~D model of a reversed fi«2d pinch. 

Here in this type of fusion device, the closed field lines 
are produced internally. A progress report on "Study of 
plasma convection and wall interactions in magnetic 
confinement system" which isaUS - department of energy 
project hrieflv explains the recently developed zero-D model. 
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appendix a 


The three equations derived in Chapter 2 can be solved 
simultaneously by the following second order method* Given a 
time step dt and the time derivation n , T , T ± at time t, then 
calculate n, T q/ at time t+dt. First calculate B, P,n Q at time 
t. Second increment n, T q/ T ± to approximate them at t+0.5 dt 
(x is any of n, T Qf 


x(t+0 *5dt) = x(t) + (0.5dt)(x c ^ 4 -^* 

t+0, 5 dt 

Third, calculate n, U g , at t+O.Sdt, using values of 
time t and n_, T , T . at t+0,5dt, Fourth, calculate T , 

G G X 6 

T . S ( I U . -T . n) /n, 
e, i 3 e, l e, i 


B, P, n Q at 
T^ from 


Finally, using n, T ,T- at t+0,5dt, increment the original 
n, T e , T i at time t to get (x is any of n, T g , T.^) 
xCt+dt) = x(t) + (dt) (x t+0 ^ 5dt ) 


APPENDIX C 

Calculation of 

Using Maxwell* s equations, 
v x B = J 

vxe= -|- 

Consider now the magnetic field diffusion equation. The plasma 
magnetic field B i (r) is related to the external magnetic 
field, plasma column area, and resistivity by a magnetic field 



diffusion equation obtained from* Maxwell's equation and 
Ohm' s law in the form 

J o = (l A i ) % + U r B ± ) 

where E is the azimuthal electric field and U is the 

r 

plasma radial velocity. Since the plasma diamagnetic current 
flows transverse to the magnetic field lines and since the 
magnetic field in the sheath is of sufficient magnitude 
so that the electron gyrofrequency is many orders of magnitude 
larger than the electron- ion collision frequency* the 
strong magnetic field limit expression for plasma resistivity 
is used in this analysis [ ^ ] 

T) = 3.27 x 10" 9 In /T?y 2 , 
j_ f e 

Reference [ ^ ] indicates that the magnetic field equation 
can be expressed as 


a V 

= vxC.t^VxB-UxB) 

3 t 

With the long compression oil length assumption, Br/Z <<1, 
the z-component can be written 


ajz 

9 t 



3 

3 r 
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r 9 r 


CrU r B z ) 


Integrating the above equation over the plasma cross 

sectional area and using the Leibnitz' s rule for differentiating 


integrals 
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Using A and U^. _ (g A/ 3t)(l/2 7i r), then the second terms 

on both sides cancel leaving 
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